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Summary 

Inviscicl aerodynamic characteristics of the Space Shuttle Orbiter were computed in support of the 
Columbia Accident Investigation. The unstructured grid software FELISA was used and computations were 
done using freestream conditions corresponding to those in the NASA Langley 20-Inch Mach 6 CF 4 tunnel 
test section. The angle of attack was held constant at 40 degrees. The baseline (undamaged) configuration 
and a large number of damaged configurations of the Orbiter were studied. Most of the computations were 
done on a half-model. However, one set of computations was done using the full-model to study the effect 
of sideslip. The differences in the aerodynamic coefficients for the damaged and the baseline configurations 
were computed. Simultaneously with the computation reported here, tests were being done on a scale model 
of the Orbiter in the 20-Inch Mach 6 CF 4 tunnel to measure the ‘deltas’. The present computations comple- 
mented the CF 4 tunnel test, and provided aerodynamic coefficients of the Orbiter as well as its components. 
Further, they also provided details of the flow field. 
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Introduction 

This report documents the results of a computational study done in support of the Space Shuttle Columbia 
accident investigation. The Columbia Accident Investigation Board report 1 concluded that the cause of the 
accident was a breach in the TPS on the leading edge of the port wing. This breach was caused by the 


1 



impact of a piece of insulation foam that had separated from the external tank soon after lift-off. During 
reentry, this breach led to the weakening of the wing, and ultimately to the loss of the vehicle. 

Examination of the Orbiter trajectory and the recorded flight data (during the re-entry until the loss of 
the signal) showed deviations in the aerodynamic coefficients as compared to the undamaged vehicle (for the 
known Mach number, angle of attack, sideslip, and control deflections). These deviations in aerodynamic 
coefficients, called the ‘deltas’, were very small at the time of Orbiter’s entry into atmosphere, but increased 
gradually reaching catastrophic values at the end. In addition, a couple of hot-spots had been noticed on 
the left side of the fuselage, possibly due to the hot gases that may have escaped from under the wing to 
the upper side. In order to investigate these observations, tests were done on a scale model at the NASA 
Langley Research Center in the 20-inch Mach 6 CF 4 Tunnel which simulates flight conditions at high Mach 
numbers. These tests were done on the baseline and several damaged configurations, and the ‘deltas’ for the 
damaged configurations were evaluated from the data. 

The purpose of the present study was to complement the CF 4 tunnel tests, and compute the ‘deltas’ 
to assist in the reconstruction of a possible damage scenario progression that would produce the observed 
variations in the ‘deltas’ with time. Since the precise damage to the Orbiter and its progression with time 
were not known during the accident investigation, several plausible scenarios that could lead to the observed 
deviations in aerodynamic coefficients were considered for the present computational study including: 

(1) Single or multiple missing RCC panel(s), 

(2) A gap behind an RCC panel that would allow hot gases to escape from wing lower surface 
to the upper surface, 

(3) Open landing gear well (due to missing landing gear door), 

(4) Dimpled wing lower surface, and 

(5) Distorted wing leading edge. 

The FELISA software was used for the aerodynamic computations, and the aerodynamic coefficients and 
the corresponding ’deltas’ were computed for the configurations listed above. One set of asymmetric flow 
computations was done to study the effect of sideslip angle on the ’deltas’. 

The FELISA Software 

The FELISA unstructured grid software is a powerful tool that can be used to build grids for complex 
geometries and compute flow solutions within a few days. This software has been used extensively in the 
past for aerodynamic computations of flow past hypersonic and planetary entry vehicles . 2,3 It has also been 
used successfully to compute high angle of attack (a greater than 40 degrees) aerodynamic coefficients of the 
Space Shuttle Orbiter . 4 The software consists of a set of computer codes for unstructured grid generation, 
and simulation of three-dimensional steady, inviscid flows using unstructured tetrahedral grids, and post- 
processing of flow solutions. Surface triangulation and discretization of the computational domain using 
tetrahedral elements are done by using two separate codes. The surface and volume grids are isotropic with 
little or no stretching. Two flow solvers are available — one applicable for transonic flows, and the other for 
hypersonic flows. The hypersonic flow solver has options for perfect gas (constant 7 ), equilibrium air, CF 4 , 
CO- 2 , and equilibrium gases representing Mars, Neptune, and Titan atmospheres. This solver also has the 
capability of solving chemical non-equilibrium flow and real gas (chemical and thermal non-equilibrium) flow. 
The hypersonic flow solver (with perfect gas air and CF 4 gas options) was used for the present computations. 
The FELISA flow solvers have been parallelized to run on multiprocessor machines. More information on 
FELISA software may be found in Peiro et al . 5 A description of the hypersonic flow solver may be found in 
Bibb et al . 6 


Geometry 

Since the Space Shuttle Orbiter has a plane of symmetry, computations were done on one half of the 
Orbiter. The ‘damaged’ Orbiter configurations, however, are all asymmetric. An examination of flow fields 
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of the baseline and damaged configurations suggested that the flow on the damaged side did not noticeably 
affect the flow on the undamaged side. The reason for this is that the damaged areas are small, and are 
located at a distance from the centerline. Therefore, it was decided to do the computations on only one half 
of the Orbiter. This not only reduced the required computer resources but also reduced the clock time by a 
factor of nearly 2. The present computations were done for the the following configurations: 

1. Baseline: A sketch of the Orbiter baseline geometry is shown in Figure 1. In order to avoid computing 
the separated flow behind the fuselage, the OMS pod, the boclyflap, and the fin, pseudo surfaces were built 
to extend these surfaces all the way to the exit plane. This causes the flow leaving fuselage and other surfaces 
not to separate, but flow smoothly over the fake surfaces to the exit plane. In hypersonic flows, this does 
not affect the upstream flow and computed ’deltas’. 

2. Missing RCC Panel(s): Figure 2 shows the ‘Panel 6 Out’, ‘Panels 6 & 7 Out’, ‘Panels 5, 6, & 7 Out’, 
and ‘Panel 9 Out’ configurations. At the time when this work was started, the actual size of the RCC panels 
and their locations on the wing leading edge were not known. Approximate location and size of panel 6 were 
extracted from an old document, and were used for the first set of computations. This is the reason for the 
panel 6 being much smaller than it should be in Figure 2a. Precise location and size of all the panels were 
later obtained, and were used for all subsequent computations. 

3. Slots: The effect of flow going through a slot behind panel 9 was studied. A sketch of the slot is shown 
in Figure 3. The length of the slot is 27.9 in. Slots of two different widths (0.53 in and 0.99 in) were studied. 

4. Open Landing Gear Well: The next damaged configuration considered was the open landing gear well 
shown in Figure 4. The well was simplified to be a rectangular hole 151 in long and 62 in wide. Three cases 
were run with 6 in, 15 in, and 30 in for the depths of the landing gear well. 

5. Dimple: The next damaged configuration studied was the ‘dimple’ (or a depression) on the underside 
of the wing near the RCC panel 9. The shape of the dimple was rectangular with the ends rounded as shown 
in Figure 5. Four dimple configurations were studied; all the dimples were 23 in wide. The first one of these 
is called the ‘Long Dimple’, and is 149.8 in long (Fig. 5a); the second is called the ‘Short Dimple’ and is 
97.8 in long (Fig. 5b). Both of these were approximately 6 in deep. The third, called the ‘Short Shallow 
Dimple’, is 97.8 in long and only 3 in deep (Fig. 5c) The fourth configuration is the ‘Long Dimple’ with the 
Panel 9 Out (Fig. 5d). 

6. Damaged Wing Leading Edge: In this case a part of the wing leading edge was assumed to be damaged 
as shown in Figure 6. The damage amounted to a change in the local angle of attack of the wing in the area 
shown. The wing leading edge was assumed to remain smooth. 

7. Full Vehicle: A set of computations were done for the full vehicle to compute the effect of asymmetric 
flow (non-zero (3). The only damaged scenario studied was the ‘Panel 9 Out’ case. 

The right-handed axis system used here has the origin at the nose of the vehicle (see Figure 1). The 
x-axis and the z-axis are in the symmetry plane; the z-axis is along the body axis and points upstream, 
and the x-axis points up. The y-axis is perpendicular to the symmetry plane and points to the starboard 
side of the Orbiter. This axis system is used for all the configurations. The reference quantities used for 
non-dimensionalizing the aerodynamic loads are as follows: 


Reference area, ( S re f ) 

Reference length for CM Y , ( c re f ) 

Reference length for CM X and CM Z , ( h re f ) 

Moment reference point 

(that is 38.5 in above the centerline, on the symmetry 


193,680 sq.in. 

474.72 in 
936.684 in 

(38.5 in, 0, -841.7 in) 
plane, and 841.7 in behind the nose.) 


Grids Used in the Present Study 

The starting point for building grids was a Space Shuttle Orbiter IGES file. The RCC panels and the 
damaged regions were defined by simple curves, and were also available in the IGES format. These IGES files 
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were processed using the software GridTool,' and a set of FELISA data files was generated. These included: 
(1) the FELISA data file that contains all the information on the surface definition and the computational 
domain in a format suitable for the FELISA grid generator, (2) the FELISA background file that specifies 
the grid spacing, and (3) the FELISA boundary conditions file. The computational domain (see Figure 7) 
was chosen to be sufficiently large so that the inflow boundaries did not affect the flow on the body. 

It was known at the outset that the ‘deltas’, which are the differences in the aerodynamic coefficients for 
the damaged vehicle and the corresponding coefficients for the undamaged or baseline vehicle were small. 
In order to capture these small quantities, large grids would be needed. Further, the grid spacing for the 
baseline and the damaged vehicles needed to be the same except in the region of the damage, so that two 
grids would be essentially the same except around the damaged region. This would insure that the resulting 
‘deltas’ would be essentially due to the vehicle damage only. 

The grid spacing was chosen to obtain dense gridding in regions of high gradients and sparse gridding 
in other places. In the freestream and in regions of very small gradients the grid spacing was 40 in, and in 
other places around the body it was much smaller (1 to 5 in). Around the fuselage nose the spacing was 
0.2 in. This provided adequate grid resolution to capture the bow shock in front of the body. Around the 
wing leading edge the grid spacing varied from 1 to 2 in. For the ‘Slot’ configuration meshes the minimum 
grid spacing within and around the slot was 0.15 in. 

Starting with the FELISA data files, unstructured surface triangulation and tetrahedral grids were gen- 
erated using the surface and volume grid generators, respectively on an SGI Onyx. A new grid was built 
for each configuration. The grids for the full configuration were especially very large, consisting of over 34 
million tetrahedral elements. The properties of all these grids are shown in Table 1. 

The grid Baseline for the baseline configuration and Baseline -6 for the panel 6 out configuration (See 
Table 1) have the same grid spacings, except around the panel 6. These grids were used for computations 
with perfect gas air. The grid Base-cff for baseline and the grids for all panel out, slot, and open landing 
gear well cases have the same mesh spacing in the computational domain except near the damaged areas. 
For the computations with the ‘Dimple’ cases, a new grid Base was built for the baseline configuration. This 
grid and all the grids used for the dimple cases (namely dimple, sdimple, ssdimple, p9dimple) and damage3 
have the same mesh spacing except near the damaged areas. 

Flow Solution 

The flow solutions were computed on a Pentium based cluster. The grids were partitioned so that the 
solution could be run on 16 or 24 processors. The FELISA hypersonic flow solver with the perfect gas air 
(7=1.4) and the CF 4 gas option were used for these computations. Each FELISA solution was started with 
solver’s low-order option, and after 2000 iterations, with solver’s higher-order option was turned on, and 
the solution was run to convergence. After every 100 iterations, the surface pressures were integrated, and 
the aerodynamic coefficients were computed. The flow solution was assumed to be converged when these 
integrated loads remained essentially constant. This required about 8,000 iterations and 200 hours of CPU 
time. 

The flow conditions used for computations with the CF 4 gas corresponded to test conditions in the 
Langley CF4 tunnel, and are as follows: 

Velocity, (Too) 2994 ft/s (=912.6 m/s) 

Density, (poo) 2.926E-05 slugs/ ft 3 (=1.509E-02 kg/m 3 ) 

Temperature, (T x ) 387.4°R (=215.2 K) 

All the computations were done for an angle of attack of 40 degrees. For the full model, computations 
were done for sideslip angles of -1, 0, and +1 degrees. 

Aerodynamic loads obtained by integrating the surface pressures were non-dimensionalized in the con- 
ventional manner, and the aerodynamic coefficients, namely Cx, Cy, Cz, CMx, CMy, and CMz were 
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obtained. The ‘deltas’ were computed by taking the differences between the coefficients for the damaged 
configurations and the corresponding coefficients for the baseline configuration. Thus, for C x , we have 

SC X 

( C x ^ damaged {C' baseline 

and similar expressions for Cy, Cz, CMx , CMy , and CMz- 

Results and Discussion 

The damaged configurations considered in the present study are such that they lead to locally large 
changes in pressure distribution on the wing leading edge. Therefore, the resulting ‘deltas’ (which are the 
difference between the aerodynamic coefficients for the damaged vehicle and the undamaged or baseline 
vehicle) are primarily due to pressure, and although the present computations are inviscid, the computed 
deltas are expected to be accurate. Further, it should be noted that these ‘deltas’ are very small. For 
example, for the Panel out cases the SCMz and the SCM x are of the order of 0.001. This numbers are 
comparable to the accuracy associated with the CMz and the CMx ■ However, as noted earlier, since the 
grids were designed specifically to capture the ‘deltas’, they are expected to be quite accurate. No attempt 
was made to made to check whether the results were grid converged. 

Out of the six ‘deltas’, the SCMx and SCMz were considered to be of primary importance. Hence, these 
two ‘deltas’ are examined in the following paragraphs. 

The Half model Results: 

Aerodynamic coefficients from the present computations for the baseline, and the damaged configurations 
on the half model are summarized in Tables 2 and 3, and the ‘deltas’ are summarized in Tables 4. The SCMx 
and SCMz for the all the damaged configurations are shown in Figures 8 and 9. It may be noted that the 
SCM\ for all the configurations are positive. Considering only the panel out cases (Figure 8), Panel 6 Out 
yields the least (0.09E-3) and Panel 9 Out gives the largest (1.42E-3) values for SCM X . The 5CM Z values 
for all panel out cases are all negative, and Panel 5, 6, & 7 Out case yields largest values (-8.0E-4). As more 
number of panels are removed from the wing, the SCMz increases. This is possibly due to the fact that as 
more area is removed from the wing there is a corresponding reduction in the normal force, that leads to 
the SCMz ■ The <5C'Mxappears to depend on the spanwise location of the panels taken out. Forward facing 
area at the wing leading edge (due to the missing panels) experiences a high pressure leading to a moment 
about the x-axis. Moving this to outboard stations leads to 5CM X ■ 

Contribution from the Slot and the Wider Slot to 5CM x and SCMz are small. 

The ‘deltas’ for the three Landing Gear Well cases show an interesting trend. The 8 CM x as well as the 
SCMz are positive for all the L.G. Well cases. The 15 in deep L.G. Well shows the largest values for both 
SCM x (1.73E-3) and SCMz (0.72E-3). For the 6 in and 15 in deep L.G. Wells the corresponding S are 
smaller. 

For the four dimple cases, as well as the damaged wing leading edge cases, SCM x as well as the SCMz 
are positive. In all these cases the wing underside was disturbed, causing the wing (underside) pressures to 
rise and produce a positive SCM x ■ It is interesting to note that the ‘deltas’ for dimple & panel 9 removed 
configuration can be obtained with a simple addition of the ‘deltas’ for the dimple alone and the ‘deltas’ for 
panel 9 removed. 

The ‘Full’ model Results: 

Computations were done for the full configuration to compute the effect of sideslip on the ‘deltas’. First, 
computations were done for the baseline configuration for /3=-l, 0, and +1 deg. and a=40 deg. with the CF 4 
tunnel conditions. Next, these computations were repeated for the damaged configuration with the RCC 
panel 9 removed from the port wing. 

The upper and lower surface C p distributions for /?= 0 deg. are shown in Figure 10a and 10b. It may be 
noted that while the absence of panel 9 changes the pressure distribution on the upper surface, whereas the 
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lower surface pressure distribution remains unaffected. The C v distributions on the body on the damaged 
and undamaged sides are shown in Figure 11a and lib. The flow going through the gap (due to missing 
panel 9) and directed towards the fuselage may be seen in Figure 10a. Impingement of such a flow on the 
fuselage side is the likely cause of the hot-spots recorded in the Orbiter flight data. 

The aerodynamic coefficients for these configurations are listed in Table 5, and are shown plotted in 
Figures 12 and 13. The ‘deltas’ are listed in Table 6. Over the /3 range of -1 to +1 degrees, all the 
aerodynamic coefficients vary linear with /?, and the ‘deltas’ are nearly constant. These ‘deltas’ (for panel 
9 configuration) computed on the full model are comparable to the corresponding values obtained using the 
half model. 
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Configuration 

Grid 

Name 

No. of 
Points 

No. of 
Triangles 

No. of 
Tets 

No. of 
Nodes 

Baseline 

Baseline 

227,478 

113,741 

14,397,123 

2,436,675 

Panel 6 Out 

Baseline-6 

230,956 

115,480 

14,502,499 

2,455,381 

Baseline 

Base.cf4 

358,092 

179,048 

17 million 

2,999,211 

Panel 6 Out 

Base.cf4-6 

347,526 

173,765 

17,261,782 

2,941,844 

Panels 6 & 7 Out 

Base.cf4-67 

340,870 

170,437 

17,420,293 

2,965,840 

Panels 5, 6, & 7 Out 

Base.cf4.567 

342,444 

171,224 

17,510,427 

2,981,318 

Panel 9 Out 

Base.cf4-9 

340,136 

170,070 

17,467,887 

2,973,650 

Slot 

slot-cf9 

425,958 

212,979 

14,044,231 

2,429,516 

Wider Slot 

Wslot.cf9 

471,274 

235,637 

22,037,778 

3,760,853 

L.G. Well, 6” Deep 

door6 

342,260 

171,132 

17,093,737 

2,913,683 

L.G. Well, 15” Deep 

door 15 

347,516 

173,760 

16,969,322 

2,894,275 

L.G. Well, 30” Deep 

doorSO 

354,408 

177,206 

17,631,557 

3,005,540 

Baseline 

Base 

237,984 

118,994 

14,638,278 

2,479,392 

Dimple 

dimple 

289,470 

144,737 

15,705,294 

2,669,783 

Short Dimple 

sdimple 

270,514 

135,259 

15,288,618 

2,595,904 

Short Shallow Dimple 

ssdimple 

264,954 

132,479 

15,218,857 

2,582,882 

Dimple & Panel 9 Out 

p 9 dimple 

270,224 

135,114 

13,035,429 

2,223,760 

Damaged LE 

damaged 

234,256 

117,130 

14,531,093 

2,460,779 

Baseline 

full~cf4 

605,450 

302,727 

34,605,278 

5,872,566 

Panel 9 Out 

full.cf9 

607,280 

303,642 

34,757,554 

5,898,267 


Table 1 Properties of grids used in the present computations. 
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Configuration 

Cx 

Cy 

Cz 

CM X 

C My 

CM Z 

Baseline 

1.2635 

1.5298E-01 

-4.7708E-02 

-2.9344E-02 

-5.4721E-02 

1.8623E-01 

Panel 6 Out 

1.2623 

1.5301E-01 

-4.8661E-02 

-2.9422E-02 

-5.5202E-02 

1.8579E-01 


Table 2 Aerodynamic coefficients for baseline and ‘Panel 6 Out’ configurations, Perfect Gas Air 
M=6, a=40 deg. 


Configuration 

Cx 

C Y 


CM X 

CM y 

CM Z 

Baseline 

1.2353 

1.5413E-01 

-5.1355E-02 

-3.0274E-02 

-2.7207E-02 

1.7772E-01 

Panel 6 Out 

1.2345 

1.5427E-01 

-5.3352E-02 

-3.0102E-02 

-2.8058E-02 

1.7753E-01 

Panels 6 & 7 Out 

1.2286 

1.5341E-01 

-6.0480E-02 

-2.9077E-02 

-2.9308E-02 

1.7633E-01 

Panels 5, 6, & 7 Out 

1.2263 

1.5363E-01 

-6.3677E-02 

-2.8128E-02 

-3.0189E-02 

1.7612E-01 

Panel 9 Out 

1.2288 

1.5475E-01 

-5.9921E-02 

-2.7435E-02 

-2.7769E-02 

1.7709E-01 

Slot 

1.2371 

1.5452E-01 

-5.1685E-02 

-2.9734E-02 

-2.9760E-02 

1.7734E-01 

Wider Slot 

1.2372 

1.5465E-01 

-5.1795E-02 

-2.9655E-02 

-2.9619E-02 

1.7735E-01 

L.G. Well, 6” Deep 

1.2468 

1.5237E-01 

-6.4568E-02 

-2.7925E-02 

-3.4559E-02 

1.7890E-01 

L.G. Well, 15” Deep 

1.2488 

1.4863E-01 

-7.4180E-02 

-2.6806E-02 

-3.4120E-02 

1.7915E-01 

L.G. Well, 30” Deep 

1.2400 

1.5371E-01 

-5.8086E-02 

-2.8897E-02 

-2.9841E-02 

1.7778E-01 

Baseline 

1.2372 

1.5425E-01 

-5.1508E-02 

-2.9758E-02 

-2.9786E-02 

1.7739E-01 

Dimple 

1.2424 

1.5255E-01 

-5.5533E-02 

-2.8742E-02 

-3.2511E-02 

1.7875E-01 

Short Dimple 

1.2398 

1.5251E-01 

-5.5451E-02 

-2.8850E-02 

-3.0767E-02 

1.7807E-01 

Short Shallow Dimple 

1.2388 

1.5342E-01 

-5.3689E-02 

-2.9295E-02 

-3.0309E-02 

1.7783E-01 

Dimple & Panel 9 Out 

1.2364 

1.5434E-01 

-6.3901E-02 

-2.5754E-02 

-3.3010E-02 

1.7824E-01 

Damaged Wing 

1.2385 

1.5741E-01 

-5.4966E-02 

-2.8747E-02 

-3.0378E-02 

1.7739E-01 


Table 3 Aerodynamic coefficients for Baseline and several damaged configurations, CF 4 gas, M=5.85 
a=40 deg. 



Model 

Gas 

SCx 

SC Y 

6C Z 

SC Mx 

SC My 

SCM Z 

Panel 6 Out 

Air 

-1.20E-03 

0.03E-03 

-0.95E-03 

0.08E-03 

-0.48E-03 

0.44E-03 

Panel 6 Out 

cf 4 

-0.40E-03 

0.07E-03 

-1.00E-03 

0.09E-03 

-0.43E-03 

-0.10E-03 

Panels 6 & 7 Out 

cf a 

-3.35E-03 

-0.36E-03 

-4.56E-03 

0.60E-03 

-1.05E-03 

-0.70E-03 

Panels 5, 6, & 7 Out 

cf a 

-4.50E-03 

-0.25E-03 

-6.16E-03 

1.07E-03 

-1.49E-03 

-0.80E-03 

Panel 9 Out 

cf a 

-3.25E-03 

0.31E-03 

-4.28E-03 

1.42E-03 

-0.28E-04 

-0.32E-03 

Slot 

cf a 

0.90E-03 

0.20E-03 

-0.17E-03 

0.27E-03 

-1.28E-03 

-0.19E-03 

Wider Slot 

cf a 

0.95E-03 

0.26E-03 

-0.22E-03 

0.31E-03 

-1.21E-03 

-0.19E-03 

L.G. Well, 6” Deep 

cf a 

5.75E-03 

-0.88E-03 

-6.61E-03 

1.17E-03 

-3.68E-03 

0.59E-03 

L.G. Well, 15” Deep 

cf a 

6.75E-03 

-2.75E-03 

-11.4E-03 

1.73E-03 

-3.46E-03 

0.71E-03 

L.G. Well, 30” Deep 

cf a 

2.35E-03 

-0.21E-03 

-3.37E-03 

0.69E-03 

-1.32E-03 

0.03E-03 

Dimple 

cf a 

2.60E-03 

-0.85E-03 

-2.01E-03 

0.51E-03 

-1.36E-03 

0.68E-03 

Short Dimple 

cf a 

1.30E-03 

-0.87E-03 

-1.97E-03 

0.45E-03 

-0.49E-03 

0.34E-03 

Short Shallow Dimple 

cf a 

0.80E-03 

-0.42E-03 

-1.09E-03 

0.23E-03 

-0.26E-03 

0.22E-03 

Dimple & Panel 9 Out 

cf a 

-0.40E-03 

0.05E-03 

-6.20E-03 

2.00E-03 

-1.61E-03 

0.43E-03 

Damaged Wing 

cf a 

0.65E-03 

1.58E-03 

-1.73E-03 

0.51E-03 

-0.30E-03 

0.00E-00 

Table 4 

The 

‘deltas’ for the damaged 

configurations at a= 

: 40 deg. 
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Config. 

P 

deg 

Cx 

Cy 


CM X 

CM y 

CM Z 

Baseline 

+1 

1.2351 

5.4679E-03 

-5.1542E-02 

-1.2911E-03 

-2.6805E-02 

1.4204E-03 

Baseline 

0 

1.2351 

-2.0009E-05 

-5.1534E-02 

3.2192E-06 

-2.6757E-02 

1.2426E-05 

Baseline 

-1 

1.2352 

-5.4895E-03 

-5.1539E-02 

1.3039E-03 

-2.6781E-02 

-1.4011E-03 

Panel 9 Out 

+1 

1.2317 

5.9294E-03 

-5.5943E-02 

2.2551E-04 

-2.7191E-02 

1.1316E-03 

Panel 9 Out 

0 

1.2317 

5.4120E-04 

-5.5853E-02 

1.5334E-03 

-2.7122E-02 

-2.4101E-04 

Panel 9 Out 

-1 

1.2318 

-4.8427E-03 

-5.5777E-02 

2.8523E-03 

-2.7124E-02 

-1.6106E-03 


Table 5 Aerodynamic coefficients for the full model, CF 4 gas, M=5.85, a = 40 deg. 


Config. 

P 

deg 

SC : x 

SC Y 

SC Z 

SC Mx 

8CM y 

5CM z 

Full 

+1 

-3.40E-03 

0.46E-03 

-4.40E-03 

1.52E-03 

-0.39E-03 

-0.29E-03 

Full 

0 

-3.40E-03 

0.56E-03 

-4.32E-03 

1.53E-03 

-0.37E-03 

-0.25E-03 

Full 

-1 

-3.40E-03 

0.65E-03 

-4.24E-03 

1.55E-03 

-0.34E-03 

-0.21E-03 


Table 6 The ‘deltas’ for the wing Panel 9 Out configuration, CF 4 gas, M=5.85, a=40 deg. 
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Baseline Conliguralion 


Fig. 1 The ‘Baseline’ configuration. 
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Fig. 2 The RCC panels. 
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Fig. 4 The landing gear well. 
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Fig. 5 The ‘Dimple’ configurations. 
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Fig. 6 The ‘Damage’ configuration. 
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Fig. 7 The computational domain. 
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Fig. 8 5CMx and SCMz for panel out, slot, and landing gear well Configurations. 
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Fig. 9 SCMx and 5CMz for Dimple and Damaged wing configurations. 
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a) Upper surface 


b) Lower surface 


Fig. 10 Wing surface C p distributions for the full model, CF 4 gas, M=5.85, a = 40 deg ., (3=0 deg. 
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Fig. 11 Body surface C v distributions for the full model, CF 4 gas, M=5.85, a=40 deg., (3= 0 deg. 
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Fig. 12 Variation of Cx, Cz, and CMy with f3 for the Baseline and Panel 9 Out configurations, CF 4 
gas, M=5.85, a=40 deg. 
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Fig. 13 Variation of CMx , CMz, and Cy with j3 for the baseline and Panel 9 Out configurations, CF 4 
gas, M=5.85, a=40 deg. 
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